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Two-component model for the deuteron electromagnetic structure 
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We suggest a simple phenomenological parametrization for all three deuteron electromagnetic 
form factors, and show that a good fit on the available data, with a minimal number of parameters, 
can be obtained. The present description of the deuteron electromagnetic structure is based on 
two components with different radii, one corresponding to two nucleons separated by ~2 fm, and a 
standard isoscalar contribution, saturated by u) and (j) mesons, only. 

PACS numbers: 


I. INTRODUCTION 

The electromagnetic structure of hadrons is characterized by a set of electromagnetic form factors (FFs), which are 
functions of one variable, the four momentum transfer squared, Q^. 

Since the pioneering work of Hofstadter Q , a large number of experimental data have been collected about hadron 
electromagnetic FFs, especially in the space-like region. It is important to note that a relatively simple phenomeno¬ 
logical parametrizations can be found for the description of the dependence of different FFs, despite the possible 
complicated dynamics which determines the hadronic electromagnetic structure. Therefore, even without a complete 
understanding of the internal structure of the hadrons, simple analytical formulas have been suggested for the hadronic 
FFs. 

For example, in case of charged pion, the corresponding FF, F.n.(Q^), can be written as: 

F^{Q^) = il + Q^/ml)-\ ( 1 ) 

where rrip is the p-meson mass, and > 0 in the space-like region. 

Similar formulas have also been suggested for FFs of the electromagnetic transition 7 -|- 7 * ^ = 7 r°, 77 , 77 ' 

and rjc'. 

Pp^^.(g 2 )=Fp( 0 )(l + 7Q"Mp)-', ( 2 ) 

where mp is a fitting parameter, which depends on the type of pseudoscalar meson. 

For a long time, the nucleon electromagnetic FFs have been described by a very simple form: 

Gpp(Q2) = GMp{Q^)/ttp = GMn{Q^)lttu = [1 + Q' [GeV2]/0.71]-2, GpniQ^) = 0, (3) 

where yip = 2.79(/i„ = —1.91) is the magnetic moment of proton (neutron). But the last experiments at the Jefferson 
laboratory, using the polarization transfer method 0 , showed a large deviation from parametrization and the 
proton electromagnetic FFs have, instead, the following behavior 

ttpGEp{Q^)/GMpiQ^) = 1 - 0.13(g2 [GeV]" - 0.04). (4) 

The experimental data about Gp„((5^), obtained from quasielastic scattering of longitudinally polarized electrons 
W a polarized deuteron target, d{e,en)p and from the measurement of the neutron polarization in d{e,e'n)p 
P, show, indeed, that GpriiQ^) 7 ^ 0, with the following parametrization |^: GpriiQ^) = —T\GMn/(X + 5. 677 ), with 
77 = ( 3 ^/( 4777 ,^), and m is the nucleon mass. 

Note that the analysis 0 of all existing experimental data concerning elastic electron-deuteron (ed) scattering 
- the structure functions A{Q^) and B{Q^) and the deuteron tensor polarization - in framework of the impulse 
approximation, leads to comparable values of Gp„(( 3 ^) and Gpp(( 3 ^) at relatively large momentum transfer, > 2 
GeV^. 
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Next in the list of simple parametrizations of hadronic electromagnetic FFs, we can mention that the magnetic FF 
for the transition N —> A(1232), = 3/2+ can also be parametrized in a simple form |^: 

Gna{Q^) = Gna{ 0){1 + QVm*^)-\ m*^ < 0.71 GeV^. (5) 

Finally, the nucleon axial FF, for the transition W* +p n {W* is the virtual IF-boson), can be described 

by the following simple formula [9j: 

G^(Q 2 ) = G^( 0 )(l + QVmi)-” ( 6 ) 

with niA = 1-06 GeV, if n = 2. 

These simple parametrizations are very useful for the discussion of different processes of elastic and inelastic scatter¬ 
ing of electrons on nuclei at high energies. Nevertheless, these parametrizations can not be considered a complete and 
precise description of the hadronic electromagnetic structure in the full region of momentum transfer (space-like and 
time-like region as well). An evident example is that the parametrizations m and 0 suggested for the description 
of the nucleon structure in the space-like region, violate the relation GsriiQ^) = GMn(Q^), at = —4 to^, i.e., at 
the threshold of the annihilation process e'^ + e~ ^ N + N. 

To avoid this and other problems, in Ref. M (updated in Ref. ^3)) another parametrization of nucleonic FFs 
has been suggested, which can be extended in the whole region of momentum transfer squared. The basic idea of 
this parametrization is the presence of two components in the nucleon structure, with different radii: the intrinsic 
structure, with radius ~ 0.34 fm (updated 0.49 fm), characterized by a dipole FF (which is the same for electric, 
magnetic, proton and neutron FFs) and a meson cloud, (which contains the p, w, (j) contributions) different for each 
of the four nucleon FFs. 

In this paper we generalize this two-component picture to the case of deuteron electromagnetic FFs, with the aim 
to find a simple parametrization for all three FFs, at least in the region where they have been completely determined, 
for < 2 GeV^. This parametrization should be useful in corresponding calculations of deuteron electromagnetic 
processes, such as e~ + d ^ e~ + d, e~ + d ^ e~ + n + p, e~ + d ^ e~ + d + etc. 

Any specific additional assumption about the validity of impulse approximation, the role and size of meson exchange 
currents, relativistic corrections, model (relativistic or non relativistic) of deuteron structure etc. is not needed. 
Instead, we parametrize the three deuteron FFs in a simple form, with a small number of parameters, normalized for 
= 0 to the electric charge, the magnetic moment and the electric quadrupole moment of the deuteron. In order to 
decrease the number of independent parameters, we will use the experimental constraint of the position of the node 
for the electric and magnetic FFs. 


II. FORMALISM 

The matrix element for ed elastic scattering, for the one-photon approximation, Fig. ^ can be written as: 

M = ■^u{k2)-ff,u{ki)J^{d), (7) 

where ki and k 2 are the four momenta of the initial and final electron, and J^^id) is the electromagnetic current of 
the deuteron. Applying the conservation of this current, the P and C invariance of the electromagnetic interaction of 
hadrons, we can write J^i{d) in the following general form 0: 

Md) =- {Gi(Q2)([/*.[/,)(p,+p2 )^+ G2{Q^)[U,^{U^-q)-U;^{U^-q)]- (8) 

G3(Q")^(C/i • g)(C/ 2 * • q){p, +P2)m} , 

where pi(Ui) and P 2 {U 2 ) are the four momenta (polarization four vectors) of the initial and final deuteron, and the 
polarization four vectors satisfy the condition: pi ■ Ui = p2 ■ U2 = 0, M is the deuteron mass. 

The deuteron invariant form factors Gi, 1 = 1 — 3, are related to the charge, quadrupole and magnetic FFs by: 

G, = Gi + ^tG2, Gm = G2, G, = Gi - G2 + (1 + t)G3, t = (9) 

The differential cross section for elastic electron deuteron scattering can be expressed in terms of two structure 
functions, A{Q^) and B{Q^), which depend on the three electromagnetic FFs: 

A(Q2) -t R(Q^)tan2 


da 

dQ 


= ctm 


( 10 ) 
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where ctm = cos^(0/2)/[4£’^ sin'*(0/2)] is the Mott cross section. Here E and E' are the incident and scattered 

electron energies, 6 is the electron scattering angle, = iEE' s\t? {9/2) is the four-momentum transfer squared and 
a is the fine structure constant, a = e^/ 47 r = 1/137. The elastic electric and magnetic structure functions and 

B{Q'^) are written in terms of the charge, quadrupole and magnetic FFs Gc{Q^), Gq{Q'^), and Gm{Q^) as: 

Am = Gim + Ir^Gim + ^TGim. sm = ^r(i+ mim- (n) 

In order to determine the three FFs, one needs another observables, usually the component t 2 o of the tensor po¬ 
larization of the recoil deuteron, in an unpolarized collision, which contains the following combination of the three 
FFs: 


~ |2''‘G'cGg + + 3 ''' [^ + ^(1 + t) tan^(0/2)] G^l , 

where S = A + B tan^(0/2). The existing data on the differential cross section m and ^20 M for electron deuteron 
elastic scatte ring allow the extraction of the three electromagnetic deuteron FFs up to ~ 2 GeV^. This has been 
done in Ref. where the world data were collected and three different analytical parametrizations were suggested, 
with a number of parameters varying from 12 to 33. 

In general, the existing parametrizations of deuteron electromagnetic FFs contain a large number of parameters, 
and are often based on analytical formulas, with poor physical content. An attempt to find a global description based 
on the vector dominance model, and satisfying the asymptotic conditions predicted by QCD at large O ^. Ilfl . lead to 
a twelve parameters fit. The fit was updated including the world data on ed elastic scattering in Ref. |1 !tI| where two 
other fits were suggested. One is a sum of inverse polynomial terms, where the first node of the corresponding FFs 
was introduced in a global multiplicative term. The number of free parameters, necessary to obtain x^/nd/ = 1.5, was 
eighteen. The last parametrization is a sum of gaussians, with some physical constraints on the parameters, which 
are the width and the position of the maximum of the gaussians. In total the parametrization contains 33 parameters 
for /ndf = 1.5. 

We suggest here a simpler parametrization, based on transparent physical content, with a minimal number of 
parameters. More precisely, we extend a two-component parametrization, already successfully applied to nucleon 
electromagnetic FFs ^3^3 recently to strange nucleon FFs [HIS , to the deuteron electromagnetic FFs. 

The deuteron is an isoscalar particle, therefore, considering only the contribution of the isoscalar vector mesons, lo 
and (/, one can write: 


Gi{Q^) = N^gi{Q^)Ei{Q'^), i = c,q,m 


( 12 ) 


with: 


— 1 — ai — Pi + at 


to: 


P^- 


to2 +Q2 + 


where rriuj {m^) is the mass of the uj ((/)-meson. Note that the dependence of Ei{Q^) is parametrized in such form 
that Fi{0) = 1 , for any values of the free parameters ai and Pi, which are real numbers. 

The terms gpQ"^) are written as functions of two parameters, also real, 7 ^ and 5i, generally different for each FF: 


g^m = ^/[^ + l^Q^X 

and Ni is the normalization of the i-th FF at = 0: 

Nc = GM = 1 , 


(13) 


Nq = Gq{0) = M^Qd = 25.83, 


M 

Nm = Gm{ 0 ) = —dd = 1-714, 

TO 


where Qd, and g,d are the quadrupole and the magnetic moments of the deuteron. 
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III. RESULTS AND DISCUSSION 


In Ref. the existing data on ed elastic scattering, differential cross section and the polarization observables, 
were reconsidered. Values of the three deuteron FFs were extracted, for the values where ^20 rneasurements were 
available, getting the values of A{Q^) and B{Q^) from an interpolation of the data on the differential cross section. 
While the magnetic FF, Gm{Q'^), is directly related to the extraction of the charge and quadrupole FFs 

requires the solution of two quadratic equations, which may lead, in some cases, to two possible roots. Therefore, the 
analysis of Gc{Q^) and Gq{Q^) consists in two different sets of solutions and two corresponding fits. 

The experimental data for Gc and Gm show the existence of a zero, for ~ 0.7 GeV^ and Qq^ — 2 GeV^. The 
constraint of a node gives the following relation between the parameters Ui and (3i, i = c and m: 


^ ml + Qh _ r. ml + Ql^ 

Ql. ‘^^ml + Qld 


(14) 


In the fitting procedure this relation allows to obtain a better description of the data and a faster convergency, 
reducing the number of free parameters. 

The expression m contains four parameters, Oi, f3i, 7 i, 6i, generally different for different FFs. We consider the 
region <2 GeV^, where the separation of Gc and Gq has been done. 

In principle the parameters 7 and S, Eq. m, may be fixed by the asymptotic behavior of the deuteron FFs, which 
follows from quark counting rules 0 . However, the range of applicability of the present parametrization is a priori 
restricted to < 2 GeV^ and this region is expected to be far from the asymptotic region Ml. 

The data basis of the present study consists in the data tabulated as in Table I of Ref. ^ and completed by 
more recent measurements from Ref. [^ . 

As a result of the procedure for the extraction of the values of GdQ^) and Gq{Qd from A(Q^), B{Q‘^) and t 2 o(Q^)j 
some experimental points show a large asymmetry of the errors, which can not be neglected in this analysis. While 
there is still no general guide how to treat asymmetric errors, we used two different ways to handle them. At first no 
asymmetry was assumed and the average of the upper ((T+) and lower (ct-) errors was taken. Then, two approaches 
recently proposed in Ref. 0 were applied. Model 1, as proposed in the paper, assumes linear dependencies and 
defines the contribution to a modified as 


X 


2 



for Ci > 0, ^ for < 0, 

i 


(15) 


where is the discrepancy between the i-th experimental point and the value of the corresponding function. Model 2, 
preferred by the author, is based on a quadratic approximation. It was not always suitable for our analysis, because, 
in some cases it doesn’t give real solutions for the contribution to 

In any case, the analysis which takes into account the asymmetry of errors (CSJ gives significant reduction of in 
all cases, but it didn’t influence significantly the resulting parameters of the fit, except for Gq(Cp), where the errors 
on the parameters were significantly reduced. 

The results were firstly obtained with a three parameter fit /3, 7 , 5, and the constraint m for Gc and Gm and a 
four parameter fit a, /3, 7 , 6, according to Eq. itT^ . for Gq. 

The parameters S and 7 are similar for all FFs Gc and Gm, with good accuracy (Table HJ. A change within 10% 
in the position of the node, slightly affects the quality of the fit, improving in general one parametrization, while the 
other gets worse. The fit is quite sensitive to the choice of initial parameters, in particular for Gq. In case of Gq, 
which is not constrained by a node, a good fit can be obtained with a large cancellation of the terms driven by a and 
/?• 

The parameters 7 and S characterize the global structure of the deuteron, and the factor 1/(1 + is related 

to the two nucleon core of the deuteron. But the isoscalar structure of the electromagnetic structure of the deuteron 
FFs, which is described by the functions Fi{Q^), in terms of (p- and cj-mesons contributions, is different for the three 
FFs, with different sets of parameters a and /3. 

From Table^one can see that the parameters 7 and 5 are not so different for the three FFs. This means that FFs 
would be mostly sensitive to the meson cloud. In order to test this, a global fit was performed, keeping the 7 and <5 
the same for the three FFs, and fitting (or fixing) a and /3 as previously. In such fit, two solutions appear also for 
Gm, related to the choice of the other two FFs. 


^ The value for Gq corresponding to Q = 2.788 fm ^ should be 2.59/g'®(, instead of 2.59(±0.073) |2^ 
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a 


7 [GeV ]-2 

5 

X^lndf 

Gc (I) 

5.9 

± 

0.1 

-5.2 ± 

0.2 

13.9 ± 1.4 

0.96 

± 

0.07 

0.8 

Gc (II) 

5.0 

± 

0.2 

-4.5 ± 

0.3 

11.5 ± 1.2 

1.11 

± 

0.09 

1.2 

G.(I) 

3.1 

± 

1.1 

- 2.1 ± 

1.2 

7.2 ±2.8 

1.6 

± 

0.5 

0.5 

G,(II) 

1.4 

± 

2.0 

- 0.1 ± 

2.4 

7.7 ± 1.6 

1.7 

± 

0.4 

0.8 

Gm 

3.78 

± 

0.04 

-2.87 ± 

0.04 

11.4 ±0.5 

1.07 

± 

0.03 

1.5 


TABLE I: Parameters for the three deuteron electromagnetic FFs. In case of Gc and Gm, ot is not a parameter, but it is 
derived from Eq. iinj. 



a 

P 

Gc (I) 
Gc (II) 

5.75 ± 0.07 

5.50 ±0.06 

-5.11 ±0.09 

-4.78 ± 0.08 

G.(I) 

G,(II) 

4.21 ±0.05 

4.08 ± 0.07 

-3.41 ±0.07 

-3.25 ±0.09 

Gmil) 

Gmill) 

3.77 ± 0.04 

3.74 ± 0.04 

-2.86 ±0.05 

-2.83 ±0.05 


TABLE II: Parameters a and /3 obtained from a global ht of the three deuteron electromagnetic FFs. The parameters 5 and 
7 are the same for all form factors and in case of Gc and Gm, a is not a parameter, but it is derived from Eq. CU- 


In Figs. mu and 21 the data points used in the present ht are shown, together with the result of this last ht 
(solid (dashed) lines correspond respectively to the hrst (second) parametrization). Open symbols in Figs. Eland 2] 
correspond to the second solution for Gc and Gq. The values of the best ht parameters are reported in Tabled The 
common parameters are 5 = 1.04±0.03, 7 = 12.1±0.5, for the hrst solution, corresponding to x2/n(i/ = 1.1, whereas, 
for the second one, 6 = 1.05 ± 0.03, 7 = 12.1 ± 0.5 and x2/nc(/' = 1.5. 

In Ref. [ 23 , in order to study the behavior of deuteron FFs, in framework of QCD, a reduced deuteron FF, fniQ^), 
was dehned as: 


= + ( 16 ) 

where TOq = 0.28 GeV^ and Fn is a generalized nucleon electromagnetic FF^. 

Comparing Eqs. m and one can identify gi(Q^) with the term (1 + Q^lml) then one expects <5=1 and 
7 = llml- From Table 1 it appears that d ~ 1, but 7 is larger and corresponds to ml ~ 0.1 GeV^. The product 
= 6jS is related to the radius of the two-nucleon (2N) component, and one hnds: i? 2 JV=l -7 fm, for TOq ~ 0.08 GeV^. 

Note, in this respect, that the standard nonrelativistic description of the deuteron results in i? ~ l/^/lmED — 4 
fm, where E]j is the deuteron binding energy. 

Data are not expected to be extended at higher in next future. So, the present parametrization can not be 
constrained at higher Q^. Nevertheless, we compared the predictions for the structure function A{Q'^), which has been 
measured up to 6 GeV^, as well as for the observables B{Q'^) and t 20 (Fig- 0- The description is good, as expected, 
in the range constrained by the fit. The difference between the experiment and the suggested parametrization visible 
for A{Q^) at larger depends on the parametrization used, due in particular to the position of a second node which 
appears for Gc{Q^)- The individual charge, magnetic and quadrupole contributions to A{Q^) are shown in Fig. 21 for 
the two possible solutions. 

The discrepancy between the data for A(Q^) and the suggested deuteron FFs parametrizations can be interpreted 
as an indication that the region at > 2 GeV^ is a transition region from the hadronic description of the deuteron 
structure to the quark degrees of freedom. 

Another possible source of the difhculty to extend the fit at higher may originate from the fact that at some 
point, at large the one-photon approximation, that is at the basis of the relations used here among FFs and 
experimental observables, and usually assumed in electron hadron elastic and inelastic scattering, does not hold 


^ A discussion of the dependence of fniQ^) on the choice of Fjv can be found in Ref. d- 
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anymore. A mechanism, where two photons, which equally share the momentum transfer squared, could become 
important. A discussion of this problem, concerning precisely ed elastic scattering data, can be found in Ref. p^. 


IV. CONCLUSIONS 

We suggested a simple parametrization of the three deuteron electromagnetic FFs, with a minimal number of 
parameters, based on a transparent physical picture. It can be used in the comparison of different theoretical models 
with experiments involving deuterons, and for a precise analytical interpolation of the experimental points in the region 
< 2 GeV^. The present parametrization is based on a classical description of the deuteron structure, in terms of 
hadronic (nucleon+meson) degrees of freedom. The region of Q^, where the separation of the charge and quadrupole 
deuteron FFs has been done, can not be easily extended in next future [^. For > 1.8 GeV^, polarization 
measurements (which normally require a secondary scattering or a polarized target) are extremely challenging at the 
present accelerators, with the present techniques of polarized targets or polarimeters due to the steep decreasing of 
the cross section. 

The parametrization proposed here can be considered a generalization of the model M , developped for the nucleon 
electromagnatic FFs. As in the nucleon case, the considered parametrization obeys, by construction, to the analyticity 
properties of FFs and can be extended to the time-like region. This is the object of a forthcoming work. 
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e^(ki) e^(k2) 



FIG. 1: Feynman diagram for electron-deuteron elastic scattering, within the one-photon mechanism. 
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FIG. 4: Fit to deuteron quadrupole form factor data. Notations as in Fig. 
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FIG. 6: Different contributions to the structure function A, according to Parametrization I (thin lines) and II (thick lines). 
The term related to Gc, Gq, and Gm, are shown as dashed, dotted, and dash-dotted lines, repectively. 












